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The measurement and analysis of electron scattering 
from ^He and ^He by Sick and collaborators reported 
20 years ago remains a matter of current interest. By 
unfolding the measured free-proton charge distribution, 
they deduced a depression in the central point nucleon 
density, which is not found in few-body calculations 
based on realistic potentials. We find that using wave 
functions from such calculations we can obtain good fits 
to the He charge distributions under the assumption that 
the proton charge size expands toward the center of the 
nucleus. The relationship to 6-quark Chromo-Dielectric 
Model calculations, is discussed. The expansion is larger 
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than than the predictions of mean field bag calculations 
by others or our CDM calculations in the independent 
pair approximation. There is interest here in the search 
for a "smoking gun" signal of quark substructure. 

I. INTRODUCTION 

The charge distribution of nuclei has been the subject of exper- 
imental studies for more than forty years. Electron scattering and 
muonic atoms now provide detailed descriptions of the full range 
of stable, and many unstable nuclides. Unique among the nuclides 
are the isotopes ^He and "^He because they exhibit a central den- 
sity about twice that of any other nuclide. There is a long-standing 
apparent discrepancy between the experimentally extracted charge 
distributions and detailed theoretical structure calculations which 
include only nucleon degrees of freedom. 

McCarthy, Sick and Whitney |l|J^ performed electron scatter- 
ing experiments on these isotopes up to momentum transfers of 
4.5 fm~^ yielding a spatial resolution of 0.3 fm. They extract a 
"model independent" charge distribution, which means that their 
analysis of the data is not based upon any assumed functional form 
for the charge distributions. Their charge distributions are shown 
in Fig. 1. Taken alone, they do not appear to be extraordinary. 
However, using a finite proton form factor, which fits the exper- 
imentally measured rms radius of about 0.83 fm, they unfolded 
the proton structure from the charge distributions to obtain the 
proton point distributions. For both isotopes there is a significant 
central depression of about 30% extending to about 0.8 fm. Sick 
also presented results where relativistic and meson effects are 
included. These are shown for ^He in Fig. 2. One note of caution 
here is that it is not possible to subtract these effects from the 
experimental data in a completely model-independent way. 
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FIG. 1. Model-independent charge distributions extracted from ex- 
periment. Reproduced from McCarthy et al.[l\. 
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FIG. 2. Point-proton density distributions obtained by unfolding the 
finite free proton form factor, allowing for meson exchange corrections 
and relativistic effects. Reproduced from Sick[2]. 

One might assume that such a central depression is to be ex- 
pected because of the short-range repulsion of the nucleon-nucleon 
interaction. So far, this is not borne out by numerous detailed the- 
oretical calculations, none of which finds a significant central de- 
pression, certainly not of the above magnitude. Relatively smaller 
central depressions are found in Green's function Monte Carlo 
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(GFMC) calculations of the alpha particle for realistic models of 
the two- and three-nucleon interaction, (see Fig. 4 below.) 

The status of theoretical structure calculations through mass 
number 4 is very satisfactory at present. Given any assumed in- 
teraction, the few body problem can be solved to within tenths of 
an MeV in energy and the wave function can be calculated to a 
precision better than that required for the present discussion. 

In using a nuclear wave function to construct a charge distribu- 
tion, one must assume a nucleon charge density and the possibility 
of meson exchange contributions. While the meson exchange con- 
tributions in the transverse channel are well-constrained (at least 
at moderate momentum transfer) by current conservation, no such 
constraint is available in the longitudinal channel. Indeed, meson 
exchange current contributions are of relativistic order and hence 
one must be careful when interpreting them with non-relativistic 
wave functions. 

Given these caveats, it is possible to reproduce reasonably well 
the longitudinal form factors of 3- and 4-body nuclei within a nu- 
cleons plus meson-exchange model. ||^J^ The current and charge 
operators are constructed from the N-N interaction and required to 
satisfy current conservation at non-relativistic order. The result- 
ing meson-nucleon form factors are quite hard, essentially point- 
like. This raises the possibility of explaining the form factors 
in quark or soliton based models, which would describe the short- 
range two-body structure of the nucleons in a more direct way than 
is available through meson exchange current models. For example, 
see the discussion of a model by Kisslinger et al. |^ below. 

We present here a possible explanation of the electric form fac- 
tors which does not involve a hole in the point-proton density dis- 
tribution, but rather is consistent with theoretical few body cal- 
culations. It involves the variation of the proton charge size as a 
function of density, or as a function of nucleon-nucleon separation. 
This is not depicted as an average 'swelling' of the nucleon, but 
as a result of short-range dynamics in the nucleon-nucleon system. 
The results presented here are preliminary but encouraging. 
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II. QUARK SUBSTRUCTURE OF NUCLEI AND 
NUCLEONS 



Within the context of sohton models, there have been numerous 
calculations of the nucleon size in nuclear media. Most of these in- 
volve immersion of solitons in a uniform (mean) field generated 
by other nucleons. Another approach has been the 3-quark/6- 
quark/9-quark bag models, which has been applied to various nu- 
clear properties, including the EMC effect. It has been applied 
by Kisslinger et al. ^ to the He electric form factors with some 



In a series of papers, Koepf, Pepin, Stancu and Wilets 
have studied the 6-quark substructure of the two-nucleon problem, 
and in particular obtain the variation of the quark wave functions 
with inter-nucleon separation. Contrary to previous expectations, 
the united 6-quark cluster does not exhibit a significant decrease in 
the quark momentum distribution function in spite of an increase 
in the volume available to the individual quarks. This is due 
to configuration mixing of higher quark states. Such a momentum 
decrease was proffered as an explanation of the EMC effect. How- 
ever, the united cluster does have approximately twice the volume 
of confinement of each 3-quark cluster, and the quarks extend to 
a volume nearly three times that of the 3-quark clusters, again 
enhanced by configuration mixing of excited states. 

In Fig. 3 we exhibit the proton rms charge radius Vp extracted 
as follows from the calculations of Pepin et al. p[: the abscissa 
gives the effective nucleon-nucleon separation vnn obtained by the 
Fujiwara transformation; the soliton-quark structure is a 6-quark 
deformed composite. The proton rms radius is defined to be 



where the quark density used in calculating < >= / pqT^d^r is 
the six-quark density normalized to unity. 



success. 
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FIG. 3. Proton rms charge radius as a function of inter-nucleon sepa- 
ration. The dashed hne is a gaussian approximation, normahzed to the 
free value. 

For well separated solitons, the tnn is just the separation of 
the soliton centers and = 0.83 fm as indicated by the horizontal 
line. Large deformations (near separation) are difficult to calculate 
so that the figure does not reproduce well the separation region. 
Shown also in the figure is a gaussian approximation fitted to rjv at = 
0, = 1 fui, and the asymptotic region, see Eq. (5) below, with 
A = 0.45, s = 0.92 fm. Then the charge distribution due to two- 
body correlations is 

fpi^i, rj; r) = [dip exp [ - |r - ri|7&^(nj)] 
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+ exp [ - |r - rj\yb\nj])}/2n'/H'(rj) 



(2) 



where we assumed and indicate explicitly that 6 is a function of the 
distance between the nucleons i and j, as we expand upon later, 
and that the proton and neutron functions are the same. Here 
"p" stands for "proton" and the Kronecker deltas pick out protons 

among i and j. 

Using the independent pair approximation (IPA) and Eq. (2) 
we find the charge distribution by employing a two-body correlation 
function p2(ri, Vj), 



Pch{r) = / d^ri / d^rj pal^i, rj)fp{ri, rf, r)/3 . (3) 



There are six pairs (i, j). Each proton appears three times; Hence 
the factor 1/3. 

III. PHENOMENOLOGICAL DENSITY-DEPENDENT 
PROTON FORM FACTOR 

To obtain some quahtative feehng for the expansion of the pro- 
ton charge size with nuclcon density, we assume a proton form 
factor fp (differing from Eq. (2)) with a size that depends simply 
on the local density and hence on the distance from the center of 
the nucleus (r') . Then 






where p is the (theoretical) point proton density 



fp{r,r') = exp[- |r 



(5) 



and we choose 



6(r')=6o[l + Ae-^"/^'], 



(6) 
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with bo = y2/3 0.83 fm, the free proton vahie. A fairly good fit to 
the data was found with A = 0.45 and s = 0.65 fm corresponding 
to 6(0) equal to the central value given in Fig. 3. The best fit, 
with only slightly better x-squared, was obtained with A = 2.10 
and s = 0.13 fm, which does not seem to be reasonable, in that the 
A is too large and the s too small. 
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FIG. 4. ^He density distributions including a fit with a variable pro- 
ton charge size. The dashed curve labelled "Carlson" is based on a 
Green's function Monte Carlo calculation [4] 



IV. THE INDEPENDENT PAIR APPROXIMATION 

In the spirit of the independent pair approximation, the charge 
distribution was calculated using Eq. (3) with p2 the two-particle 
correlation function The proton size parameter 6(ri2) was first 
taken from the gaussian fit to the calculations of Pepin et al. 
The improvement over the free constant proton size, as shown in 
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Fig. 5 (dot-dash) was small. 

A phenomenological fit to the data was made with a parame- 
terized b given by 



6(ri2) = 6o[l + ^e-('--/^)"], 



(7) 



where n = 2 is of the gaussian form, n > 2 yields a sharper 
transition. Indeed, n — oo yields a step function. Recall that the 
model of Kisslinger et al. corresponds to a step function. 

We examined n — 2, A and 6. Although the A and s were dif- 
ferent in each case, the quality of fits were very similar. The corre- 
sponding best fit values of {A, s, n) for the three n's were (2.185, 
0.883, 2), (0.976, 1.245, 4), (0.774, 1.34, 6). In Fig. 5 we show the 
results for n=2 since the others are indistinguishable to the eye. 



0.20 



0.15 — 




a 




point density, Carlson p 
free proton 

Pepin 
2.185, .883, 2 
charge density, expt 



0.05 



0.00 



FIG. 5. *^He density distributions constructed from variational point 
densities and two-body correlation functions in a parameterized varia- 
tional calculation by Carlson et al.[4]. The curve labelled "Pepin" uses 
the gaussian fit of Fig. 3, based on the calculations of Pepin et al.[8]. 
The solid curve is a phenomenological fit as described in the text 
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V. CONCLUSIONS 



We have obtained a phenomenological fit to the "^He charge dis- 
tribution by assuming a proton size which increases with increasing 
density. More specifically, we minimize the charge distribution 
squared using a two-parameter gaussian function of r, the distance 
from the center-of mass. 

We would like to identify the variable proton size with the struc- 
ture function of Fig. 3 derived from 6-quark A^-A^ studies in the 
spirit of the independent pair approximation. Fairly good agree- 
ment with experiment was obtained with a phenomenological pa- 
rameterization of the proton size function. 

The inadequacy of the previous calculation might be due to 

• A constant confinement volume was assumed for the six-quark 
structure as a function of deformation. It may be that the inter- 
mediate volume (between separated and united clusters) is larger. 

• The independent pair approximation may be invalid at the 
high densities of the central region. 

• Meson effects should be recalculated using the quark structure 
functions given (say) by the six-quark IPA model. 

Items 1 and 3 are topics for further investigation. In addition, 
one must study the predictions of such models for quasi-elastic scat- 
tering. In the quasi-free regime, nucleon models produce a good 
description of the data as long as realistic nucleon interactions, in- 
cluding charge exchange, are incorporated in the final-state interac- 
tions. Unlike the charge form factor, two-body charge operators 
are expected to play a much smaller role here [4]. The combination 
of the two regimes provides a critical test for models of structure 
and dynamics in light nuclei. 
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